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0\ \ We present a model calculation of the generalized parton distributions where the nu- 
cleon is described by a quark core surrounded by a mesonic cloud. In the one-meson 
approximation, we expand the Fock state of the physical nucleon in a series involving a 
^ . bare nucleon and two-particle, meson-baryon, states. We discuss the role of the different 
Fock-state components of the nucleon by deriving a convolution formalism for the unpo- 
CN I larized generalized parton distributions, and showing predictions at different kinematics. 

O 

\0 • 1. The meson-cloud model for the nucleon 
O 

r-| ■ The convolution model for the physical nucleon, where the bare nucleon is dressed by its 
virtual meson cloud, has a long and successful history in explaining properties such as form 
Qh! factors [H] and parton distributions In this paper it has been revisited and applied 
^ I for the first time to study Generalized Parton Distributions (GPDs) that have recently 
been introduced and discussed in connection with Deeply Virtual Compton Scattering 
(DVCS) and hard exclusive meson production (for reviews, see Refs. [ El IH ISl E] ) • 
The basic assumption of the meson-cloud model is that the state of the physical nucleon 
5^ I N can be decomposed according to the meson-baryon Fock-state expansion as a super- 
position of a bare nucleon state and states containing virtual mesons associated with 
recoiling baryons. This state, with four-momentum = Patj,) = {Pn,Pn) and 

helicity A, is an eigenstate of the light-cone Hamiltonian 

Hlc = E [Hiiq) + H^{q) + Hj{N, BM)] . (1) 

B,M 

In Eq. (P), HqIq) stands for the effective-QCD Hamiltonian which governs the constituent- 
quark dynamics, and leads to the confinement of three quarks in a baryon state; analo- 
gously, ifQ^(g) describes the quark interaction in a meson state, and Hi{N, BM) is the 
nucleon-baryon-meson interaction, and the sum is over all the possible baryon and meson 
configurations in which the nucleon can virtually fluctuate. In the calculation presented 
here we use perturbative treatment of the meson effects. So we truncate the Fock space 
expansion of the nucleon state to the one-meson components, and expand the nucleon 
wavefunction in terms of the eigenstates of the bare Hamiltonian Ho = HqIq) + HQ^{q). 
The corresponding state of the physical nucleon |A^) can be written as 

Ip., A; N) = Vz\p,, A; AT) + E / E k.) 

B,M ^l^TTj ^yi^l^y^yy, 
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X kx + ypN±, A'; 5) 1(1 - -kx + (1 - y)pN±, A"; M), (2) 

where we introduced the function (f)x!'yf'^^^\y, kj_) to define the probabihty amphtude for 
a nucleon with hehcity A to fiuctuate into a virtual BM system with the baryon having 
hehcity A', longitudinal momentum fraction y and transverse momentum k^, and the 
meson having helicity A", longitudinal momentum fraction 1—y and transverse momentum 
— k^. Furthermore, in Eq. Q, the constant Z, giving the probability that the physical 
nucleon is a bare core state, ensures the correct normalization of the nucleon wave function: 

p1. A'; H\p+, p^A; H) = 2{2'nfp+5{p'+ - p+)5^'\v'^ - v±)hw. (3) 



2. The unpolarized generalized parton distributions 

In the definition of GPDs we choose a symmetric frame of reference where the virtual 
photon momentum and the average nucleon momentum = \{p^i^ + p%) are collinear 
along the z axis and in opposite directions. Furthermore, = —q^q^ is the space-like 
virtuality that defines the scale of the process, t = = (p'^ — p^)^ is the invariant 
transferred momentum square, and the skewness ^ describes the longitudinal change of 
the nucleon momentum, 2^ = — A+/p^. 

For each flavor q the soft amplitude corresponding to unpolarized GPDs reads 



A^) = /^e^^^>-(p^,A'^|^(-|z)7+^(Hb^,A^) 



(4) 



where x defines the fraction of the quark light-cone momentum [k = xp%), Xn (A^) 
is the helicity of the initial (final) nucleon, and the quark-quark correlation function is 
integrated along the light-cone distance z~ at equal light-cone time {y~^ = 0) and zero 
transverse separation (zj_ = 0) between the quarks. The leading twist (twist-two) part of 
this amplitude can be parametrized in terms of the chiral-even helicity conserving GPD, 
H'^{x,C,, A_|_), and the helicity fiipping GPD, E'^{x,C,, A^), for partons of fiavor q. 

In the following we will derive the convolution formulas for the GPDs in the three 
different regions corresponding to ^ < a; < 1, — ^ < a; < ^, and — 1 < x < — ^. 

2.1. The region ^ < x < 1 

In this region the GPDs describe the emission of a quark from the nucleon with mo- 
mentum fraction x + ^ and its reabsorption with momentum fraction x — ^. In the 
meson-cloud model, the virtual photon can hit either the bare nucleon or one of the 
higher Fock states. As a consequence, the DVCS amplitude can be written as the sum of 
two contributions 

^AlA.(^,e, A^) = ZF0-(x,e, A^) + 5F,\^,^(x,e, A^), (5) 



where F'^''^"'^^ is the contribution from the bare proton, described in terms of Fock states 
with three valence quarks, and 6F'^ is the contribution from the BM Fock components of 
the nucleon state, corresponding to five-parton configurations. The valence-quark contri- 
bution Fy^'^^^ can be calculated in the light-front overlap representation derived in Ref. [ 
[7] , and applied to the case of iV = 3 valence quarks in in Ref [ El E] , where one can also 
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find the explicit expression in terms of bare-nucleon light-cone wave functions (LCWFs) 
derived in a constituent quark model. The 5Fy term in Eq. (0) can further be split 

into two contributions, with the active quark belonging either to the baryon (SF''^^^^) or 
to the meson i.e. 

SFl.^ix, e, A^) = E [6Fll^ (X, A^) + SF'J'[^{x, ^, A^)] . (6) 

B,M 

The first term in Eq. © corresponds to the case when the baryon is taken out from 
the initial proton with a fraction ys + ^ of the average plus- momentum pj^, and after 
the interaction with the initial and final photons is reinserted back into the final proton 
with a fraction ys ~ ^ of the average plus- momentum p^. The transverse momentum of 
the baryon is psj, — Aj_/2 before, and p^_|_ + A_l/2 after the scattering process. The 
meson substate is a spectator during the whole scattering process. As final result, the 
convolution formula for 6Fl!^J^ reads 

where Fy^ is the scattering amplitude from the active baryon in the BM component of 
the nucleon, and can be obtained in terms of the baryon LCWFs as explained in Ref. [ 

US]. 

Analogously, we can derive the meson contribution to the scattering amplitude, corre- 
sponding to the case when the pion takes part to the interaction process while the baryon 
remains as a spectator. In such a case the role of the meson and baryon substates is 
interchanged with respect to the situation described before. The meson is taken out from 
the initial proton with a fraction yj^j + ^ of the average plus-momentum p^, and after the 
interaction with the initial and final photons is reinserted back into the final proton with 
a fraction yj^.^ — ^ of the average plus- momentum p^. The transverse momentum of the 
meson is Pm± ~ A_|_/2 before, and Pm± + A_|_/2 after the scattering process. Viceversa, 
the baryon substate is inert during the whole scattering process. Therefore the meson 
contribution to the scattering amplitude is given by 



X 



B A,A',A" 

^\n{N,BM)i^ ~ ^^^X'j,{N,BM) 



1 - m, -kM±) [0^/)^; ' ^(1 - Vm. -kM±)]*, (8) 



where Fy^^ is the scattering amplitude from the active meson in the BM component of 
the nucleon, given explicitly in Ref. [E] in terms of meson LCWFs. 

2.2. The region -1 < ^ < 

In this region, the scattering amplitude describes the emission of an antiquark from the 
nucleon with momentum fraction — and its reabsorption with momentum fraction 
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— (x — ^). Therefore, the only nonvanishing contribution comes from the active antiquark 
in the meson substate of the BM Fock component of the nucleon wavefunction, i.e. 

A^) = A^), (9) 

where SFy^^{x, ^, A_|_) is given by the same convolution formula (jHI), with the integration 

range over y^j between —x and 1, and with the LCWF overlap representation of Fy^ in 
the range — 1 < x < —C,- 

2.3. The region —^<x<^ 

In this region, the scattering amplitude describes the emission of a quark-antiquark pair 
from the initial proton. In the Fock-state decomposition of the initial and final nucleons 
we have to consider only terms where the initial state has the same parton content as the 
final state plus an additional quark-antiquark pair. In the present meson-cloud model, 
the initial state is given by the five-parton component of the \BM) Fock state, while 
the final state is described by the three-valence quark configuration, multiplied by the 
normalization factor In principle, we can distinguish between two cases: i) the 

active quark belongs to the baryon substate, and the active antiquark is in the parton 
configuration of the meson substate of the initial nucleon; ii) both the active quark and 
antiquark belong to the meson substate of the initial nucleon and the baryon is a spectator 
during the scattering process. However, this last contribution is vanishing because it 
involves the overlap of two orthogonal states, i.e. the wave functions of the baryon in 
the initial state and of the bare nucleon in the final state. As a consequence, the only 
non vanishing contribution corresponds to the case i), and it can be explicitly derived in 
terms of the overlap between the five and three parton components of the LCWFs of the 
physical nucleon state as explained in details in Ref. [ITH]. 



3. Results 

In this section we present results for the GPDs in the meson cloud model by restricting 
ourselves to consider only the pion-cloud contribution and disregarding the contributions 
from mesons of higher masses which are suppressed. As a consequence, the accompanying 
baryon in the {Bn) component of the dressed proton is a nucleon or a A. For the bare- 
hadron constituents of the nucleon state we use the LCWFs in the relativistic constituent 
model adopted in previous works to describe the valence contribution to GPDs both in 
the chiral-even and chiral-odd sector [ |Sl IHl HDl HH 1121 US]- The model assumption to 
describe the vertex functions (px^^/^^^ are given in details in Ref. [I13j. 

First let us study the forward limit, ^ = 0, t = 0. In this limit the scattering ampli- 
tude without nucleon helicity flip reduces to the ordinary parton distribution, and the 
convolution formulas for H'^{x, 0, 0) coincides with the formulation of the meson cloud 
model for the ordinary parton distributions in deep inelastic processes [ |2] . In Fig. ^ the 
spin-averaged H"^ and the helicity-fiip E'^ GPDs are plotted together with the separated 
contributions from the bare proton (dashed-dotted line), the baryon (dashed lines) and 
the meson (dotted lines) in the baryon-pion fluctuation. All these contributions add up 
incoherently to give the total result (full curves). The bare proton contribution is always 
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Figure 1. The different contributions to the spin-averaged {H'^, upper panels) and helicity- 
flip {E'}, lower panels) generalized parton distributions calculated in the meson-cloud model 
for flavours u (left panels) and d (right panels), at ^ = and t = 0. Dashed lines: baryon 
contribution from the \BM) component. Dotted lines: meson contribution from the \BM) 
component. Dashed-dotted lines: contribution from the bare nucleon. Full curves: total result 
as a sum of the different contributions. 



positive within its support (0 < x < 1) with the exception of E'^ for which it is neg- 
ative. The same behaviour characterizes the baryon contribution from the baryon-pion 
fluctuation that is also limited to the range < x < 1, consistently with the assumption 
that the only active degrees of freedom for such a baryon are the valence quarks. The 
sea-quark contribution, extending all over the full range — 1 < x < 1, is determined by 
the antiquark residing in the meson of the baryon-pion fluctuation. The resulting effect of 
the pion cloud is thus to add a contribution for negative x and to increase the magnitude 
of the GPDs for positive x with respect to the case of the bare proton. In particular, for 
positive and small x the pion cloud contribution as a whole is comparable to that of the 
bare proton, confirming the important role of the sea at small x. We also note the faster 
fall off of E'^ with respect to if^ for x — > 1, showing the decreasing role of the Melosh 
transform to generate angular momentum in E"^ with increasing quark momentum. 

In all cases at x = the GPDs have a zero. This is due to fact that in the overlap 
integrals the various terms of the proton wavefunction are taken at one of their end points. 

The forward limit of the first moment sum rule for the spin-averaged GPDs, is cor- 
rectly fulfilled. For the helicity-fiip GPDs, the first moment sum rule gives the quark 
anomalous magnetic moments, for which we find the values = 1.14 and k'^ = —1.03. 
Although these numbers are 5% and 10% off the experimental values for the up and down 
quarks, respectively, we found that the contribution of the pion cloud gives a substantial 
improvement with respect to the results obtained in the model with only valence quarks. 

Going beyond the forward limit, we show in Fig. El the results for the GPDs at ^ = 
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Figure 2. The different contributions to tlie spin-averaged {H'^, upper panels) and lielicity- 
flip {E'^, lower panels) generalized parton distributions calculated in the meson-cloud model for 
flavours u (left panels) and d (right panels), at ^ = and t = —0.2 GeV^. Line style as in Fig. ^ 



at t = —0.2 GeV^. The relative contribution of the different components is not modified 
by switching on the momentum transfer t, only the overall magnitude is decreased. This 
is in agreement with the common believe that the main part of the t dependence of the 
GPDs is exhibited by their first moments, i.e. by the quark Dirac and Pauli form factors. 

In Fig.iniare plotted the isoscalar and isovector combinations of GPDs plotted at ^ = 0.1 
and t = —0.2 GeV^. We see that GPDs in the ERBL region are rather regular functions 
over the whole range, with zeros at the endpoints x = This result is quite different 
from the oscillatory behaviour predicted by the chiral quark-soliton model [ ^ where 
the valence contribution of the discrete level is a smooth function extending into the 
ERBL region and adding to the sea contribution. Here this is forbidden because the 
support of the valence contribution is limited to the DGLAP region. In addition, the 
transition amplitude between the bare-proton and the baryon-meson components vanishes 
at X = ±C,, approaching these points both from the ERBL and DGLAP region. This 
generates a discontinuity of the first derivative of GPDs at x = ±^ which, however, is not 
in contradiction with general principles [EJ- 

Furthermore, in the DGLAP region both for x > ^ and x < — ^ no striking difference 
arises in Fig.Elfor the spin-averaged GPDs H'^^'^ with respect to the results in the forward 
limit shown in Fig. ^ while for the helicity-flip GPDs the (negative) d contribution coming 
from the baryon in the \BM) component is responsible for a broader shape ai x > ^. 

4. Conclusions 

We described a convolution model for the unpolarized GPDs where the nucleon is viewed 
by a quark- valence core surrounded by a mesonic cloud. This meson-cloud model gives the 
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Figure 3. Isoscalar (u + d, left panels) and isovector (u — d, left panels) combinations of the 
spin-averaged (upper panels) and helicity-flip (lower panels) generalized parton distributions 
calculated in the meson-cloud model, at ^ = 0.1 and t = —0.2 GeV^. Line style as in Fig. ^ 



possibility to link GPDs calculated in the light-front formalism to the nucleon description 
in terms of constituent quarks including a sea contribution already at a low-energy scale. 
Since the contribution in the ERBL region is vanishing in the forward limit, it can not 
be easily inferred from parametrizations in terms of parton distributions. Therefore, the 
present calculation gives new insights to model the off-forward features of the GPDs, and 
can be further used as a suitable input at the hadronic scale to study the behaviour under 
evolution at higher scales. 
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